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1. INTRODUCTION 

Segmentation is a science in the field of digital image processing that aims to find more information 
contained in images. In the medical field, segmentation can assist medical personnel in identifying areas of 
the patient’s body with the help of an magnetic resonance imaging (MRI) or computerized tomography (CT) 
machine [1]. By utilizing the right segmentation method, the disease can be detected accurately and making it 
easier for medical personnel to take steps for healing or adventure [2]. The technology used to describe the 
inside of the human body includes MRI and CT. Each technology has advantages and disadvantages. MRI is 
safer than the CT method because it does not emit radiation, but the results of the MRI machine image are 
affected by noise from the scan results [3]. MRI images can provide a lot of useful information if this 
information can be extracted correctly and correctly. For example, in the MRI image of the brain which 
contains information about the tumor using a precise and accurate segmentation method, the area of the 
tumor can be extracted accurately as well. Brain tumor segmentation has become a popular research topic in 
recent years [4]. 

Several segmentation methods that have been proposed and conducted include the use of popular 
clustering algorithms such as the fuzzy C-means (FCM) clustering algorithms [5], [6] and the use of classical 
clustering methods such as K-means combined with other algorithms such as in [7] which combines K-means 
with mean shift segmentation to reduce file size and detect brain tumors. From all of the clustering methods 
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that have been proposed and conducted, there are still some weaknesses, namely the inability to accurately 
distinguish brain tumor areas due to inhomogeneity problems on MRI images. Therefore, the brain tumor 
segmentation method developed by utilizing the neural network method to detect tumors as in research [8], [9]. 
From the results obtained from the segmentation using the neural network, it can be concluded that the 
method proposed in each study obtained satisfactory results, but the weakness of each method that uses a 
neural network is the relatively long computation time. Therefore, a segmentation method based on the 
region/area of the image has emerged, called region-based segmentation. Region-based segmentation 
methods such as the active contour model (ACM) are one of the most widely used methods in brain tumor 
segmentation research [10], [11]. 

In general, the ACM method is divided into 2, namely edge-based and region-based, the edge-based 
segmentation model is based on information from the edges of the image to determine the path of the 
algorithm to the edges of the expected image, whereas the region-based model, segmentation runs based on 
information from the edges of the image. Information contained in the surrounding area, by calculating the 
statistical value of each region in the image. Many studies that use the ACM method have been carried out, 
including by Banday and Mir [12], who proposed a semi-automatic brain tumor segmentation method where 
user interaction is required to determine the initial area of the ACM method, uses the gray level co-ocurrence 
matrix (GLCM) and gray level run-length matrix (GLRLM) features to generate internal forces and external 
forces for ACM, and the Gabor filter is based on an edge-based ACM model. The pricipal component 
analysis (PCA) method is also used to reduce the features that have been extracted so that 3 features of each 
extraction feature are selected as parameters in the research. Another study using the ACM segmentation 
method was proposed by Pham et al. [13], who proposed 2 fitness functions to segment brain tumors 
accurately, taking into account local information from MRI images which also remains reliable in dealing 
with bias in MRI images. Sun et al. [14] also researched the segmentation of brain tumors using a 
combination of 2 algorithms. In his research, we used fuzzy-based energy calculations which were used as 
parameters for the ACM algorithm so that the curve would not easily fall to the local minima value. Another 
combination method proposed by Sheela overcomes the initial initialization problem of the ACM method and 
also the edge area problem of the segmentation results. Sheela and Suganthi [15] uses radius contraction and 
expansion to select the initial initialization area, then the FCM algorithm is used to eliminate the resulting 
areas from ACM that still have non-tumor areas. Research related to the Chan-Vese-based ACM method was 
also conducted by [16] who makes a system that can determine the most appropriate method for segmenting 
brain tumors based on previously extracted features, in the form of 12 key features including the mean, 
harmonic mean, trimmed mean, max value, median, standard deviation, skewness, kurtosis, entropy, contrast, 
energy, and homogeneity, and it was found that the Chan-Vese ACM method can be relied on in brain tumor 
segmentation. Combined methods for segmenting brain tumors are proposed by [2] which is divided into 2 
stages, first one is eliminating the part of the normal brain using the threshold function then using the Harris 
extraction feature algorithm which detects the edges and corners of the observation area to determine the 
greatest energy value. From the results of the extraction features, the convex hull method is used to roughly 
estimate the portion of the brain tumor which is then used as an initial initialization for the ACM algorithm. 
Another research that using the combination of several algorithm such as Kapur’s entropy and Chan-Vese 
algorithm conducted by provide 2 proces to segment brain tumor in MRI images. Bat algorithm and Kapur’s 
entropry are used for extracting the trilevel of threshold from the image, then segmentation process is done 
using active contour segmentation and also as post-processing to identify the boundary of the tumor area. 

The use of the brain tumor segmentation method based on the region was also developed by [17] by 
adapting the Bayesian theorem method as a parameter to determine the relation between pixels which can 
solve the inhomogeneity intensity problem that exists in the MRI image, besides that the proposed method is 
claimed to be able to segment tumors accurately even though it is affected by bad noise, by using the Markov 
random field as a function of energy with how to form a windows area which will be used to calculate the 
relation between pixels. The window area size used in this study is 8. The development of the segmentation 
method based on the active contour model is becoming more advanced, as in [3], which developed the 
extended local patch fuzzy ACM with a modified superpixel method which is expected to increase the 
robustness of the proposed method. The extended localized patch-based fuzzy active contour (ELPFAC) 
method that is proposed itself relies on the observation area which is divided into 2 parts which then 
calculates each energy from each part with additional energy calculations in the part that has a black area 
such as the image background, this is because in the previous method, which was developed by the author, 
the black area/background image that is included in the observation area makes a different representation of 
the energy in the gray part of the image, in this case, the brain tumor area, so an ELPFAC method was 
developed to solve this problem. The ACM method can indeed provide good segmentation results if the 
image to be segmented has homogeneous pixels, but in fact, this is not possible, therefore Sandhya et al. [18], 
proposes a new method, namely a combination of self organization map (SOM) and ACM where SOM has a 


TELKOMNIKA Telecommun Comput El Control, Vol. 20, No. 4, August 2022: 825-833 


TELKOMNIKA Telecommun Comput El Control Oo 827 


role to provide information on inhomogeneity intensity for the ACM algorithm. From several studies that 
have been done it can be concluded that the ACM method is a powerful method in segmenting MRI images, 
but from the overall method a few of the methods that have been proposed mention the selection of a starting 
point as the initialization of the ACM algorithm, which is the first step in the segmentation process with the 
method. Sauwen et al. [19] proposed a semi-automated method for segmenting tumor region using regularized 
non negative matric factorization which calculate the approximate the next area from the user-selected seed for 
a fine tumor segmentation, then followed by morphological post-processing to remove the outliers from the 
segmentation result. A hybrid method for segmenting brain tumor was proposed by [2], using Harris feature 
extraction and convex hull for determining the first area initialization then followed by active contour model 
to explore the tumor region and after that a morphological operations was used to remove holes or unwanted 
pixel that still persist from the segmentation process. 

In this study, we propose a method that can provide initial initialization for the ACM algorithm in 
order to obtain accurate and relatively fast tumor segmentation, requiring as fews iteration as possible in the 
ACM algorithm. By utilizing the multi-level Otsu threshold applied to MRI images with fluid attenuated 
inversion recovery (FLAIR) modality, to find the initial initialization area, the Chan-Vese ACM algorithm is 
then used to walk to the surrounding area to find specific tumor areas. The advantages offered by our 
proposed method are the speed and low computational cost in segmenting MRI images. 


2. RESEARCH METHOD 
2.1. Preprocessing 

Before the image enters the segmentation stage, the image is subjected to a preprocessing stage to 
improve image quality. In this study, we used 2 preprocessing steps, namely noise removal and contrast 
adjustment. We also perform image cropping to acquire only the meaningful image information through 
these stages: 
— Cropping image and convert to grayscale space 

The size of the MRI image used in this study is 512x512 pixels, but the dimensions of the image 
still leave a slice notation for each extracted MRI image in the lower-left corner, both for MRI images and 
ground truth images. Therefore, we cropped the image for both types of images in order to provide accurate 
and precise segmentation results, namely in the tumor area alone without any interruption from other parts. 
Here, we obtain an image dimension of 317x311 pixels for both types of images. After the image has been 
successfully cropped, the MRI image is then converted to a gray image to maximize the results of the 
multi-level threshold algorithm. 
— Noise removal using Wiener filter 

MRI images are often used as objects of research to detect diseases in the human body, but this is 
often difficult to do because MRI images are impacted by noise due to the conversion process of the MRI 
machine readings [20]. The use of a Wiener filter is considered to be able to eliminate noise contained in 
MRI images. Here, wiener filter worked but still maintain the quality at the edges of the image, which plays 
an important role in knowing the boundary between the tumor area and the normal brain area. 
— Contrast adjustment 

The use of the multi-level threshold algorithm works well if the contrast in each image area has a 
different intensity value because this will affect the results of the multi-level threshold algorithm. Therefore, 
to maximize this, we apply a contrast adjustment algorithm that increases the contrast value of the image to 
1% at the lower limit and 1% at the upper limit of the image intensity value. This makes light areas of the 
image, such as the tumor area, more visible visually, and dark areas such as the background image darker. 


2.2. Multi-level Otsu thresholding 

The method for selecting the initial initialization area for the ACM algorithm uses the multi-level 
thresholding Otsu proposed by [21]. The multi-level thresholding method is one of the most basic methods 
for image segmentation, by utilizing the optimization of the Otsu thresholding criterion calculation [22], [23]. 
In this study, the number of threshold levels is set at 3 levels and applied to the MRI image that has been 
converted into a grayscale image. Multi-level thresholding is based on calculating the minimum value of 
variance and the maximum value of variance of an image, the initial step is to divide the image histogram 
into 2 parts, then calculate the variance value of the 2 classes to obtain the maximum value of the two classes, 
and that is what is used as the threshold value. For the multilevel Otsu thresholding formulation, an image A 
has a gray intensity value of {0,1,2,3,..., L-1}, as if a threshold value ¢ divides the image into foreground (C,) 
and background (Co), and for each member of the gray intensity values of each class are Cy = [0,1,2,3,..., L-1] 
and C,= [t,t + 1,t+2,t+3,..., L-1]. Given p; as a probability representation of the gray intensity value i that 
appears on each image pixel. The equation for calculating the probability of the background as in (1) and 
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foreground as in (2) image A [23]. To calculate the mean of the foreground and background the following as 
in (3) and as in (4); and for the foreground and background calculation equations are as shown in (5). 


0 = Dice pi (1) 

01 =i pi (2) 
Mo (t) = Diz git (3) 
m (O = Died 2 (4) 
o? (t) = wo(t).@4(t). (u(t) - w(t)” (5) 


The equation for multi-level threshold is almost similar to the Otsu threshold equation, but with the 
addition of the variable s, which contains several levels of variables t = t4, tz, t3,...,¢s, which divides the 
image into s+1 areas Cy, C,, C2, ...,Cs,,. The addition of the s variable changes the equation from 
calculating the variance between class s + 1 to be as follows as in (6). And to get the maximum t value 
which is the threshold value is as follows [21] as in (7). 

Hs = Dede bo (6) 


i= ts—1 wj 


(ti, t3, t3, t= vette ett o? (ty, tz, tz, ts) (7) 

Based on Figure 1, the red box being the non-tumor area and the blue box being the tumor area. 
We can conclude this because we have observed all the datasets we use and from all image datasets we can 
make the decision that the largest area that has been detected is a strong candidate as the initial seed for the 
initialization of the ACM algorithm. We apply morphological area filtering to select areas from the results of 
converting the MRI image to a binary image using a threshold level to extract the area which will be used as 
the ACM initialization. 


2.3. Active contour model 

In this study, the active contour model used is the Chan-Vese model developed by Chan in his 
research active contour without edges [24]. Segmented has a severe noise disturbance and inhomogeneity 
intensity in the MRI image pixels [18], [25]. This causes the segmentation results to get poor performance 
because the Chan-Vese (CV) algorithm is difficult to find the minimum energy point used to generate 
internal forces that make the contour moving outward run continuously to find the minimum energy. 
However, it does not cover up that ACM is a good method of segmenting MRI images if it is applied to the 
right image type. Some of the advantages between ACM and other segmentation methods according to [12] 
robust against noise, robust against false edges, and information on shape and texture will improve 
segmentation results [18], [26], [27]. For the CV-ACM algorithm mechanism itself is to minimize the energy 
function [2] as in (8), where C is a curve that goes along iteration, inside (C) is the inner area of curve C and 
outside (C) is the outer area of curve C, and i; also iz are input constants for the inside and outside curves, 
respectively. Chan-Vese optimized as in (9). 


F (C) + F(C) = Sinsiaeql 0) - Paxdy + J, ssiaec ll 9) - fp Padxdy (8) 


F (ix, ip, C) = a. Length(C) + p. Area(inside(C)) + kı Lelie) =i \Pdxdy + 
kz Joutsiaeco CP) — i, | dxdy m 


It is intended that the calculation of the C curve function can be calculated at the object boundary, 
where a = 0,6 = 0,k, and k, = 0, which are the constant values. The next step for the segmentation stage 
is to find the tumor area on the MRI image using the ACM algorithm. At the outset of the explanation of this 
study, we stated that our method has the advantage of being the time it takes to achieve the most optimal 
results. For this, we conducted several experiments by changing the maximum number of iterations required 
for the ACM algorithm to be able to fully segment brain tumors, and we found that the number of iterations 
of 20 times is the final result we get from several experiments that have been conducted. 
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In Figure 1, are the results of morphological area selection at the multi-level thresholding stage and 
then added morphological enhancements that have been used in research [20] which aims to improve image 
quality in the sense of connecting between pixels on a separate threshold image and eliminating pixels that 
are not included in the structural element used. As shown in Figure 1, the blue square used for initial area of 
ACM, but there is also a red square that indicate other candidate for the initial area of ACM, the parameter 
we used for choosing the optimal candidate is the biggest blob with the highest intensity value. 
The parameters in this study use a ‘disk’ shape with a radius of 5 for structural elements which aims to 
connect pixels and a ‘disk’ shape with a radius of 3 for structural elements which aims to eliminate pixels 
that are not within the radius of the structural element, so that the initialization area obtained becomes 
smoother than the threshold image. It does not affect the process of segmenting the brain tumor area. 


Multi threshold 


O m o 


Figure 1. Example results using multi-level threshold and initial area of ACM 


3. RESULTS AND ANALYSIS 

Our research uses MRI images obtained from the multimodal brain tumor image segmentation 
benchmark (BRATS) 2015 image dataset which is widely use in the brain tumor segmentation such as 
in [2], [19], [28]. We extracted 86 MRI images from 86 different patients with FLAIR modality as shown in 
Figure 2(a) along with ground truth images with slices that correspond to one another as in Figure 2(b). 
The choice of FLAIR modality MRI image data is because the FLAIR type has brighter characteristics of the 
tumor area compared to other areas, such as normal brain areas and skull bones. The dataset of MRI images 
obtained from BRATS 2015 has undergone preprocessing elimination of skull bones so that no additional 
preprocessing steps are required for this. 


(a) (b) 


Figure 2. Brain tumor segmentation: (a) MRI Image and (b) ground truth image 


The evaluation used in this study used the dice similarity (DS), Jaccard index (JI), and Hausdorff 
distance (HS) calculation which compares dice similarity between original image as shown in Figure 2(a) 
MRI image and Figure 2(b) ground thruth image with results of the segmentation image. DS and JI are for 
measure overlapping areas of the two images being compared while HS calculated the maximum distance 
between two images. The image resulting from segmentation (S) will be compared with the ground truth 
image (G) as shown in Figure 2(b) using (12). The more the DS and JI value approach 1 while HS value 
approaches 0, it can be interpreted that the image of the segmentation results is similar to the ground truth 
image, which means that the segmentation gets accurate results. 
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In addition to measuring the accuracy of segmentation, in Table 1, we also measured the time 
required for an MRI image to be segmented. All experiments were carried out using a laptop with an AMD 
FX 2.58 Ghz processor with 16 GB RAM, and Matlab 2017b. We present the data for several samples of 
MRI image data in the following Table 2. It can be seen from the results of time observations in the Table 1, 
the time needed to segment 1 MRI image is 1.5 seconds + 0.3. For the DS results, we obtained satisfactory 
results with the average DS for the entire image of 0.7905, an average of 0.6765 for the Jaccard index, and an 
average of 33.3947 for the Hausdorff dist. This can be achieved because the maximum number of iterations 
is only 20 times. As explained in the sub-chapter of the proposed method, the multilevel threshold algorithm 
changes the MRI grayscale image to a binary image with a certain threshold value, in this case, there are 3 
levels and what we use is level 3, so only part of the MRI image has a grayish intensity value more than the 
threshold level 3 will be converted into a binary image obtained binary image results which will then be 
selected to be the initial initialization point for the ACM algorithm. The binary image selection result of the 
conversion using the threshold value has the role of eliminating small pixel clusters around the largest area so 
that finally an optimal initial initialization area is obtained. In this study, we also conducted several tests for 
the maximum number of iterations required and we obtained the following data in the Table 3. 


Isnc| 


Dice similarity = amen (12) 
2 
; _ |snG} 
Jaccard index = A (13) 
duausaory (G.S) = max {max minila — DII, max minllb — all} (14) 


Table 1. Result of segmentation using ACM 


Original image Segmented image 


Name Intial mask Ground truth 


MRI (1).jpg 


MRI (2).jpg 


MRI (9).jpg 


MRI (15).jpg 


MRI (16).jpg 


Please note that we did not try all sample points for the maximum number of iterations, we used an 
interval of 10 starting from 10 to 100. From the results we obtained, it was concluded that the number of 
iterations = 20 obtained the highest dice similarity and Jaccard index value among all experiments performed. 
Mean while for HS value the best value is obtained at iterations = 50. The MaxIter obtained can be reduced 
because the process of initializing the area/starting point of the ACM algorithm is close to the actual tumor 
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area, this can be seen in the experimental results Table 1. Chan-Vese’s own ACM method is also reliable in 


terms of segmenting the tumor area as in Table 4. At the multilevel threshold stage, the addition of a contrast 
adjustment process also plays an important role in determining the area/initialization point. 


Table 2. Time taken, DS value, JI value, and HS value 


Name Time taken (s) Dice similarity — Jaccard index __ Hausdorff dist 
MRI (1).jpg 1.488436 0.7760 0.6694 13.0384 
MRI (2).jpg 1.756735 0.9097 0.8523 11.4018 
MRI (9).jpg 1.427263 0.9151 0.8344 10.0499 
MRI (15).jpg 1.556853 0.9449 0.8925 28.2843 
MRI (16).jpg 1.462690 0.9156 0.8392 15 
MRI (21).jpg 1.774595 0.9633 0.9291 10.2956 
MRI (30). jpg 1.866347 0.9075 0.8306 6.3246 


Table 3. Number of iterations with best result 
Num of MaxIter Dice similarity Jaccard index Hausdorff dist Time taken (s) 


100 0.753289 0.6328116 37.244438 99.634963 
90 0.7598404 0.6402436 35.873814 92.432717 
80 0.7661297 0.6475711 34.317097 90.242293 
70 0.7721321 0.6546112 33.388309 80.650162 
60 0.7776351 0.6607835 32.713718 73.543678 
50 0.7827799 0.6665037 32.379356 58.926755 
40 0.7870497 0.6712934 32.462101 57.151485 
30 0.790093 0.6752079 32.856724 47.900759 
20 0.7905458 0.6765056 33.394741 36.461008 
10 0.7828362 0.6680661 34.716537 33.653083 


We also compare the performance of our algorithm if we do not include the additional contrast 
adjustment stage, and we get the results for DS = 0.7166 with the time taken is = 40.206 at MaxlIter = 20. 
This proves that the addition of the contrast adjustment stage has an important role in determining the 
area/initialization point, with the difference in DS = 0.0739 and the difference in time taken = 3.745 seconds, 
which is quite a big difference. For comparison, we compared the performance of our method with several 
methods that have been proposed in Table 4. Please note that we did not try all sample points for the 
maximum number of iterations, we used an interval of 10 starting from 10 to 100. 

Based on Table 4, iterations = 20 obtains the highest dice similarity and Jaccard index value among 
all experiments performed. The MaxIter obtained can be reduced because the process of initializing the 
area/starting point of the ACM algorithm is close to the actual tumor area, this can be seen in the 
experimental results Table 4. Chan-Vese’s own ACM method is also reliable in terms of segmenting the 
tumor area. 


Table 4. Segmented image comparison 


Segmented image Ground truth Dice similarity 
e 0.64 


0.75 


Authors 
Pei et al. [28] 


Shivhare et al. [29] 


Proposed 0.79 
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4. CONCLUSION 

This study focuses on fast and accurate segmentation of brain tumors. The MRI image used in this 
study was obtained from the BRATS 2015 dataset by extracting 86 MRI images that had tumors and ground 
truth images for each patient. Our proposed method focuses on selecting the initial initialization area/point 
for the ACM algorithm. The initial initialization area selection uses the multi-level thresholding algorithm 
that implements Otsu thresholding. For the threshold level used in this study, we chose 3 threshold levels, 
with the level we used was level 3. After that, the MRI image was converted into a binary image with the 
selected threshold and then extracted the largest area as the initial initialization area/point. The next step is 
segmentation using the ACM algorithm with the parameter MaxlIter = 20. The results obtained from the 
method we propose are measured by dice similarity of = 0.7905, Jaccard index of = 0.6765 and Hausdorff 
dist of = 33.394, with the time taken = 36.461 seconds. From what we have done in this research, we can 
conclude that our proposed method is more accurate and faster than other method we compared with. 
For further research, we will develop an ACM algorithm that is more accurate in segmenting MRI images, 
while still minimizing the time taken algorithm, so that an accurate and fast, and reliable method of brain 
tumor segmentation is obtained. 
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